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ABSTRACT 

Observations with the IRS spectrograph onboard Spitzer have found many sources with very deep Si features at 9.7 /im, 
that have optical depths of r > 1. Since it is believed that a few of these systems in the local Universe are associated 
with Compton-thick active galactic nuclei (hereafter AGN), we set out to investigate whether the presence of a strong Si 
absorption feature is a good indicator of a heavily obscured AGN. We compile X-ray spectroscopic observations available 
in the literature on the optically-thick (r9.7 Mm > 1) sources from the 12 /^m IRAS Seyfert sample. We find that the 
majority of the high-r optically confirmed Seyferts (six out of nine) in the 12 fim sample are probably Compton-thick. 
Thus, we provide direct evidence of a connection between mid-IR optically-thick galaxies and Compton-thick AGN, 
with the success rate being close to 70% in the local Universe. This is at least comparable to, if not better than, other 
rates obtained with photometric information in the mid to far-IR, or even mid-IR to X-rays. However, this technique 
cannot provide complete Compton-thick AGN samples, i.e., there are many Compton-thick AGN that do not display 
significant Si absorption, with the most notable example being NGC 1068. After assessing the validity of the high 
9.7 frni optical-depth technique in the local Universe, we attempt to construct a sample of candidate Compton-thick 
AGN at higher redshifts. We compile a sample of seven high-r Spitzer sources in the Great Observatories Origins Deep 
Survey (GOODS) and five in the Spitzer First-Look Survey. All these have been selected to have no PAH features 
(EWe.2/jm <0.3/im ) to maximise the probability that they are bona-fide AGN. Six out of the seven GOODS sources 
have been detected in X-rays, while for the five FLS sources only X-ray flux upper limits are available. The high X-ray 
luminosities (Lx > 10 42 ergs -1 ) of the detected GOODS sources corroborates that these are AGN. For FLS, ancillary 
optical spectroscopy reveals hidden nuclei in two more sources. SED fitting can support the presence of an AGN in 
the vast majority of sources. Owing to the limited photon statistics, we cannot derive useful constraints from X-ray 
spectroscopy on whether these sources are Compton-thick. However, the low Lx/Le^m luminosity ratios, suggest that 
at least four out of the six detected sources in GOODS may be associated with Compton-thick AGN. 

Key words. X-rays: general; X-rays: diffuse emission; X-rays: galaxies; Infrared: galaxies 
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1. Introduction 

Hard X-rays (2-10 keV) are extremely efficient in detect- 
ing AGN, as they can penetrate large columns of dust and 
gas. This allowed the Chandra and XMM-Newton missions 
to map the AGN universe with unprecedented detail. In 
particula r, about 90% of the X-ray background h as been 



resolved (Alexander et al. 



2003 



Luo et al. 2008), reveal 



ing a sky density in the CDF-N of > 5000 deg z ( Bauer et 



al 



2004). The majority of these sources are obscured AGN 
presenting column densit ies > 10 22 cm -2 (e 



2006] |Akylas et al. 2006 [) 

Hi 



Tozzi et al 



owever, even the hard X-ray surveys may be miss- 



sources, the Compton-thick AGN, which have column den- 
sities > 10 24 cm~ 2 . Alth ough this population remains elu- 



sive (see Comastri 2004 for a review) , there is concrete ev- 
idence for its pres ence. The peak of the X-ray background 



at 20-30 keV 



Frontera et al. 2007 Churazov et al. 



2007 Moretti et al.l|2009|) can be reproduced only by invok 



mg a significant number of Compton-thick sources at mod- 
erate redshifts. However, the exact density of Compton- 
thick sources required by X-ra y background synthesis mod- 
els st i ll remains an open issue ([Gilli, Comastri k. Hasinger . 
20071 iSazonov et al.[ |2008l ITreister, Urry Virani[ |2009p. 



ing a substantial fract ion of the most heavily obscured h]&± holeg in the bcal Umverse (see [gjfe 



Additional evidence of a numerous Compton-thick popu- 
lation comes from the directly measu red space den sity of 

" 19821. 



It 



Send offprint requests to: I. Georgantopoulos, 
ioannis . georgantopoulosOoabo . inaf . it 



e-mail: found that the black hole space density is a factor of 1.5-2 
higher than that predicted from the X-ray luminosity func- 
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tion dMarconi et al. 2004] [Merloni & Heinz} 2008}. The 

the 



exact number depends on the assumed efficiency of the con- 
version of gravitational energy to radiation. Direct searches 
in the ultra-hard 20-70 keV band by Swift and INTEGRAL 
did not detect large numbers of Compton-thick sources 



(c 


g. Ajello et al. 2008 Tueller et al. 




2008 


Paltani et 


al. 




2008 Winter et al.| 


2009). Nevert 


leless, it is possi- 



5Ie that even these ultra-hard surveys are biased against 
the most heav ily obscured reflectio n-dominated Compton- 
thick sources ( |Burlon et al. 20101. Owing to the limited 
imaging capabilities of these missions, the flux limit probed 
is very bright (~ 10~ n erg cm~ 2 s _1 ) allowing only the 
detection of AGN in the local universe. At higher red- 
shifts, some Compton-thick AGN have been reported in 
the deepest XMM-Newton and Chandra observations in the 



Ch andra deep fields flTozzi et al.[|2006[|G eorgantopoulos et 
aT] |20"09j |Comastri et al.[|2011[ |Feruglio et al.[|201ip . 



Mid-IR wavelength observations have attracted much 
attention because they provide an alternative way of detect- 
ing heavily obscured systems. This is because the absorbed 
radiation by circumnuclear dust is re-emitted in the IR part 



of the spectrum. [Martmez-Sansigre et al. ( 2005 ) argue that 
a population of bright 24 /im AGN with no 3.6 /im detec- 
tions is as numerous as unobscured QSOs at high redshift 



(z > 2 ) . On the basi s of X- r ay stacking analysis, Daddi et al. 
P007|);|Fiore et al.| |2008|); |Georgantopouios et al. 



Treister et al.| ( |2009b| ); |Eckart et al.| ( |2010p ; |Donley et al. 



(2010) propose that a fraction of infrared excess, 24 (im 



bright sources are associated with Compton-thick AGN. 
These sources are found at high redshift (z ~ 2), and their 
contribution to the X-ray backg round is expected to be 
small (<1%; [Treister et al.]|2009| . In contrast, the bulk of 
the contribution to t he X-ray backgro und is produced at 
redshifts z ~ 0.7- 1 ( |Gilli et &\\ [2007] . 

At high redshifts, mid-IR spectroscopy with Spitzer- 
IRS, has detected a number of sources with large columns 
of obscuring material, as inferred from their 9.7 Si fea- 



tures (r > 1; hereafter called high-r sources; Dasyra et al 
20091. These systems could be associated with Compton 



thick AGN . A local analog sy stem is the nearby ULIRG 
NGC6240 ( |Armus et aT||2006| ), which is well kno wn in X 



ray w a veleng ths to host a Compton-thick AGN ( Vignati 
et al. 1999). There is little information available at X- 



ray w avelengt hs for these high-r systems at higher redshift 



(Bauer et al. 2010) 



The primary goal of this paper is to investigate whether 
an efficient way to identify Compton-thick AGN is indeed 
to look for sources with deep silicate absorption at 9.7 /im. 
The structure of the paper is as follows: we first compile 
a local s ample of high opt ical-depth (r g .7 Mm > 1) sources 
from the |Wu et all (|2009h flux-limited Spitzer-IRS obser- 



vations of the 12 /zm Seyfert sample of |Rush et al. (1993) 



We investigate the X-ray spectra available in the literature, 
to deduce how many of these high-r objects are heavily 
obscured or Compton-thick. Finally, we compile a sample 
of high-r AGN at higher redshift with available X-ray ob- 
servations, using the Chandra Deep Fields, as well as the 
Spitzer FLS sample. Although these sources are faint in 
X-ray wavelengths (and thus it is difficult to derive with 
great certainty their X-ray spectral properties), we attempt 
to estimate whether they are heavily absorbed from their 
L x /L6 pm luminosity ratio. 

We adopt H D = 70kms~ 1 Mpc" 1 , tt M = 0.3, and Q A = 
0.7 throughout the paper. 



2. Sample selection 

2.1. The local sample 

We employ a sample of local galaxies with a deep Si absorp- 
tion feature in the mid-IR, that are known to definitely host 
an AGN and also have X-ray data, so that the fraction of 
Compton-thick AGN among them can be quantified. An 
ideal data set fo r this purpose is the Spitzer IRS spectro- 
scopic sample of Wu et al. ( 2009 1 . These authors present 

35 / im spectra for 103 galaxies 
This is 



1993J. 

local Seyfert 



low-resolution Spitzer 5.5 
from the 12 /jm Seyfert sample (iRush et al. 
a complete, unbiased, flux-limited sample of 
galaxies, selected from the IRAS Faint Source Catalog by 
means of optical spectroscopy. 

To identify sources with high optical depth at 9.7 /im 
(rg.7) due to Si absorption, we use the 6 and 13 /im contin- 
uum flux to interpolate the intrinsic AGN unob scured flux 
at 9.7 /im in a similar way to Spoon et al. (2007). We subse- 
quently measure rg.7 by calculating the natural logarithm 
of its observed to intrinsic value. For screen extinction, this 
corresponds to —rg.7. For sources without spectral cover- 
age at 13 /im, we use the 6 and 7 /im continuum values to 
extrapolate at 9.7 /im. We reject all sources with optical 
depths smaller than one. Nonetheless, various spectral fit- 
ting techniques can lead to different rg.7 values due to dif- 
ferent continuum assumptions, see e.g. the small differences 
in o ur me asurements of rg.7 and those derived by |Bauer et| 
ah] p010) in the FLS sources. F or this reason, we also run 
TAHFTT (see [Smith et "alj |2007[ ) for all sourc es for which it 
i s poss ible, using the parameters presented in |Dasyra et al.| 
( 2009J) . We confirm that all of the sources selected with our 
primary technique also hav e rg.7 > 1 based on PAHFIT. 
The above re-analysis of the Wu et al. ([2009 ) Spitzer spec- 
tra shows that rg.7 > 1 for eleven AGN (see Table]!]). The 
IRS spectra of man y of these sources were also discussed in 
Armus et al.|(|2007|. 



2.2. High-r 9.7 /im AGN at higher redshift 

To identify sources in the distant Universe that are deeply 
obscured in the mid-IR, we used all the Spitzer-IRS spec- 
tra available for the 4 deg 2 field of the First Look Survey 
(FLS), as well as for the CDF-N and CDF-S (~900arcmin 2 
in total) of the Great Observatories Origins Deep Survey 
(GOODS). The two surveys arc complementary in target- 
ing sources for our analysis. The FLS is a shallow 24 /im 
survey, whose spectra were limited t o a depth of 0.9 mJ y 



(Yan et al. 2007| |Sajina et al.| |2007 Dasyra et al. 2009), 



targeting luminous high-z galaxies. This is ideal for iden 
tifying many high-r systems, as the highest amounts of 
obscuration are fo und in the sources with the highest in - 
frared luminosity (llmanishi, Maiolino & Nakagawa 2010). 
However, the GOODS area has the most sensitive mid-IR 
and X-ray observations available. We select again sources 
with rg.7 Mm > 1, using the methodology described above. 

Since the initial sample selection is made at infrared 
wavelengths, we expect to find a significant number of 
non-AGN sources, which owe their infrared emission to 
dust heated by star- format ion processes. To minimize 
the contamination by pure star-forming systems, we ex- 
clude sources with significant PAH emission at 6.2 or 
11.3 or /im (see Fig.fTI) , since it has been demonstrated 
that AGN lack strong PAH features (e.g. Genzel et al. 
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Table 1. Spitzer IRS properties of the high-r 9.7/im AGN in the 12/itm sample 



Name 


z 


Type 


T 


log[^L 6M mJ 


AGN 


EW(6.2) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Mrk 938 


0.020 


Sy2 


1.2 


10.09 





0.440 


NGC1125 


0.011 


Sy2 


1.0 


9.12 


0.81 


0.258 


108572+3915 


0.058 


Sy2 


3.5 


11.51 


0.98 


<0.021 


UGC5101 


0.039 


Syl 


1.4 


10.52 





0.229 


NGC 3079 


0.004 


Sy2 


1.3 


9.10 





0.458 


Mrk 266 


0.028 


Sy2 


1.0 


10.03 





0.608 


Mrk 273 


0.038 


Sy2 


1.7 


10.51 





0.192 


Arp 220 


0.018 


Sy2 


2.4 


9.86 





0.344 


1 19254-7245 


0.062 


Sy2 


1.2 


11.03 


0.88 


0.064 


NGC 7172 


0.009 


Sy2 


1.9 


9.75 


0.2 


0.045 


NGC 7582 


0.005 


Sy2 


1.0 


9.56 





0.274 



The columns are: (1) Name; (2) Redshift; (3) Optical AGN type; (4) Optical depth at 9.7pm; (5) Logarithm of vL v IR 
monochromatic luminosity; at 6/im in units of solar luminosity; (6) Fraction of AGN contribution at 6/mi according to the 
spectral decomposition based on broad-band mid-IR photometry available from NED (see sectio n 3.2.3 for details of the models 



used); (7) Equivalent-width of PAH 6.2 /im feature in units of /im taken from Wu et al. (20091 



Table 2. X-ray properties of the high-r 9.7 /mi AGN in the 12 /urn sample 



Name 


z 


N H 


log[L x j 


X-ray Rcf. 


Mission 


Comment 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Mrk 938 


0.020 


40 


42.30 


1,16 


XMM 




NGC 1125 


0.011 




41.97 


2 


SWIFT 




108572+3915 


0.058 




41.30 


3 


Chandra 




UGC5101 


0.039 


140 


41.67 


4,6,7,8,11,16 


Chandra/ XMM 


a 


NGC 3079 


0.004 


200 


40.25 


5,16 


XMM 


a 


Mrk 266 


0.028 


>160 


41.7 


1,6 


XMM/ Chandra 


a 


Mrk 273 


0.038 


40 


42.40 


10,16 


Chandra 




Arp 220 


0.018 


> 100 


40.96 


7,12,14 


XMM 


b,c 


1 19254-7245 


0.062 


> 100 


42.57 


7,8,13 


BeppoSAX 


a 


NGC 7172 


0.009 


8 


42.20 


9 


BeppoSAX 




NGC 7582 


0.005 


160 


42.00 


9,15,17 


BeppoSAX 


a 



The columns are: (1) Name. (2) Redshift. (3) X-ray column density in units of 10 22 cm 2 . (41 Logarithm o f the obs cured X-rav 

]Gv. 



luminosity i n t he 2-10 keV band in u nit s of ergs . (5) X-ra y spectroscopy reference: l)|Guainazzi et al. 1(2005 



et (|20i0[> 3)|Iwasawa et al.|fl2009|) 4) |Imanishi et al.|j2003|> 5)|Akylas fc Georgantopoulos 



7) llwasawa et al.| Q2005p 8) |Braitoet al.|(|2003|) 9) |Dadma| (120080 1U 



12 

Nandra 



Ptak ct al 



KOTTT 



;2009) 6) Brassington ct al 
11) |Gonzalez-Martin et al. 



llwasa wa et a lJ 

(120021) 15y[Bianchi et al.| (|2009|) 16) |Bngntman~ 




]2003T l3) |Georgantopoulos et aL] ( 2010p 14) ]Clements et al.| ^ 

(6) X-ray mission. (7) X-ray criterion on which the source is classified as Compton-tnick: a. detection of the 
absorption turnover b. high-equivalent width FeKa line c. Flat spectrum (see text for details). 



1998) |Hernan-Caballero et aU |2009| ). We define as PAH- 
poor sources (mostly AGJN dominated) those that have 
EW(6.2/mi)<0.3/mi, or EW(11.3/mi)<0.3/mi where the 6.2 
A wavelength is not covered. 

There are 220 sources with IRS spectroscopy in the FLS 
sample, and we end up with 20 high-r AGN in FLS us- 
ing the method described above. Among these, there are 
six sources with Chandra observations available. Five of 
them have no obvious PAH features (see Fig. [I] and Table 

t, while one (FLS-283) has EW(11.3/zm)=0.68/mi and is 
erefore not included in our sample. Pope et al. (in prepa- 
ration) have consistently reduced and analyzed all pub- 
licly available IRS spectra in both GOODS-N and ECDFS 
fields. For this paper, we search through the 150 GOODS 
IRS spectra database to find AGN with high r g .7, select- 
ing 15 sources. Eight of these sources are excluded because 
they have prominent PAH features with EW larger than 
0.3/im (see Table[5]). The properties of the 12 sources in the 
high-redshift sample (seven GOODS and five FLS sources) 



are given in Table [3j while their IRS spectra are shown in 

Fig. OD 

2.3. Comparison between the high and low redshift samples 

It is important to highlight the differences between the se- 
lection at low and high redshift. In the former, all sources 
are bona-fide AGN according to optical spectroscopy. Many 
sources would be classified as composite AGN/star- forming 
systems, on the basis of the EW of the 6.2/im PAH feature. 

Since there is no optical spectroscopy available for many 
sources in the high redshift sample, the AGN classifica- 
tion is based on the absence of strong PAH features. As 
discussed in section 4.3.2, our AGN selection criterion is 
shown to be effective for at least eight out of twelve sources 
based on their high X-ray luminosities, and optical spectra 
and most probably all sources based on the measured MIR 
excess in their SED fitting. Therefore, the local and the 
high-redshift samples have different relative contributions 
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- 0.0002 - 



» (Mm) 

FLS8245 





Fig. 1. IRS spectra of the GOODS and FLS sample. 



to the star- formation in the total IR energy budget. Four 
sources in our local sample (Mrk938 NGC3079, Mrk 266, 
Arp220) would not satisfy the EW(6.2/im) < 0.3/t/m selec- 
tion criterion imposed in the high redshift sample. 

3. Analysis 

3.1. The local sample 
3.1.1. X-ray observations 

X-ray photometric observations exist for all sources in the 
local AGN sample. For one source (NGC 1125) only Swift- 
BAT o bservations are availab le at very high energies (> 
15 keV; Cusumano et al. 2010[ ), and no X-ray spectrum can 
be derived. For another source (IR AS 08572+3915) ther e 



are Chandra observations available (Iwasawa et al 
but unfortunately they are not sensitive enough to a^ 



2009) 
low the 



For the remaining nine sources, we compile the X-ray spec- 
troscopic results available in the literature in Table[2| 

3.1.2. X-ray spectra 

The classification of a source as Compton-thick is based on 



one o f the following criteria (see Georgantopoulos et al. 
20091: 



derivation of the X-ray column through X-ray spectroscopy. 



(a) The most reliable criterion is the detection of the 
absorption turnover at high energies. These are the 
transmission-dominated Compton-thick sources with 
relatively 'mild' column densities of a couple of times 
10 24 cm -2 . In these sources, we can directly view 
the obscured component through the t orus, even 
at energies below lOk eV (e.g. NGC3079; |Akylas fc 
Georgantopoulos 2009). 
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(b) The detection of a high (~ 1 keV) equivalent-width 3.2.2. X-ray spectra 
(EW) Fe line. 

(c) A 'hard' spectrum with r ~ 1 or flatter. This is con- 
sidered to be the signature of a reflection dominated 
Compton-thick source, where the X-ray emission comes 
sole ly from reflection from the backside of the torus (see 



e.g. Murphy fe Yaqoob 2009). 



Six of the nine sources for which X-ray spectra can be de- 
rived probably host a Compton-thick nucleus based on one 
of the above criteria (see Table[2]). Out of these, UGC 5101 
is probably the most debatable case: in a study of ULIRGs 
with la rge amounts of molecular gas, |Georgantopoulos et ah] 
(2010) detect a flat spectrum (r ~ 1.2) using Chandra, to- 
geth er with a high EW (~3.6 keV) FeKa line, in agreement 
with Ptak et al. ( |2003 ). The combination of a flat spectrum 
with a high EW FeKa line suggests a Compton-thick AGN. 



Imanishi et al. (2003) analyzed both the XMM-Newton 
and the Ch andra data and found a consid erably lower EW 
(^400 eV). |Gonzalez-Martin et al.| ( |2009| ), using Chandra, 
claim the direct detection of a mildly Compton-thick AGN 
with a column density of Nh ~ 1- 4 ± 0.4 



Finally, Brightman & Nandra (2010) report a column den- 
sity of 10 z3 cm - ^ by re-analyzmg the available XMM- 

the same authors report a 
10 23 cm" 2 ) for Mrk266 



cm 

Newton data. We note that 
much lower column density (2 
(NGC5256) than the one derived by Guainazzi et al. (2005 1. 
NGC 7582 is one of the few known 'changing-look' ACN, i.e. 
sources with past X-ray observations in both a Compton - 



thin and a reflection-dominated state (Bianchi et al. 2009 ). 
The three high-r AGN that are definitely not Compton- 
thick are Mrk938, Mrk273, and NGC 7172. These have 



column densities between Nh = 8 x 10 cm and 4 



For Mrk273, Iwasawa et al. (20111 note that 



the X-ray emission may come from a different location than 
the IR nucleus. They also point out that the IR nucleus may 
be associated with a Compton-thick source, although the 
evidence for this remains weak. 



3.2. Sources at high reds hi ft 



3.2.1. The X-ray observations 



The FLS has few (and shallow; 



30 ksec) X-ra y data 
points, which are presented in Bauer et al. (20101. Their 
X-ray observations targeted primarily ULlROs. There are 
X-ray observations available (all of them yielding upper lim- 
its) for only five out of the 20 high-r AGN in FLS. However, 
the GOODS survey for the Chandra Deep Fields has the 
deepest X-r ay images to date. The CDF-N was observed 
for 2 Ms (see Alexander et al.[|2003[ ). We examine the 2 Ms 
CDF-N X-ray images and find that all sources are detected 
in the total ( 0.3-8 keV) band. Three of the sources are in 
the catalog of Alexander et al. ( 2003 1 , while the other one is 
marginally below the detection threshold of these authors. 



The CDF-S was originally observed for 2 Ms (Luo et al. 
2008), while the field was observed for an additional 2 Ms 



Here, we use the entire set of 4 Ms images. In the CDF-S, 
two sources (out of the three high-r sources in this field) 
are detected in the total (0.3-8 keV) band. 



Among the X-ray detections, we can derive low-quality X- 
ray spectra in five cases but in only three we can constrain 
both the photon-i ndex and the h ydrogen column density. 
We use XSPEC flArnaud[ |1996] fo r the spectral fitting. 
We use the C-statistic ( Cash| 19791, leaving the data un- 
grouped. The errors correspond to the 90% confidence level. 
We fit a simple power-law with photoelectric absorption 
from cold material in all cases. First, leave both the photon 
index and column density as free parameters for the three 
sources with the highest quality photon statistics. The re- 
sults are given in Table|4] There is no clear evidence of a 
Compton-thick source. GNTRS-30 presents a very flat spec- 
trum, with r = 0.54, but very large errors. To obtain more 
stringent constraints on the absorbing column density, we 
fix the p hoton-index to the commonly observed value of 
r = 1.8 ( |Dadinaj [2008| . We see that the spectra of the 
sources GN_1RS-19 and GNJRS-30 are heavily obscured 
with column densities 



3.2.3. Spectral energy distributions 

To estimate the total infrared emission of the sources in the 
high redshift sample, and assess the possibility that they 
host an AGN, we fit their broad-band spectra with com- 
binations of starburst and AGN templates. The starburst 
te mplates we use come from the SWIRE template library 
F ( Polletta et al. 2007 ) , which is a compilation of oberved 
sED o f near by galaxies, and from the sample of | Chary fe 
|Elbaz| (|2001[). The AGN tem plates we use come from Silva, 
Maiolino & Granato ( 2004 1 , who combine nuclear SEDs of 
Seyferts with a range of absorption columns. We also con- 
struct a s ample of AGN templ ates using the type-1 QSO 
SED from Richar ds et al.|fl2006 ) and a pplying dust absorp- 
tion following [Rosenthal et al* (2000), to account for the 
9.7 /im absorption feature. 



The photometry we use for the GOODS-N and GOODS- 
S sources comes from the IRAC (3.6, 4.5, 5.8, and 8.5 /im) 
public GOODS data release and the SIMPLE survey 
(gpzfe e rlRA C/MUSYC Public Legacy in ECDF-S;[Damen 



et al.| 2010) respectively. MIPS photometry at 24 and 



70 /im, comes from the FIDEL survey (Far-Infrared D eep 
Extragalactic Legacy survey; see Magnelli et al.||2009[) and 
the 85 /xm photometry in the GOODS-N from |Pope et aL 
p005l ). The FLS sources have IRAC and MIPS ( 70 /xm) 
photometric values published in Sajina et al. (2007), while 
their MIPS ( 24 /im ) and 1.2 mm fluxes are published in 



pajina et al. (2008). In the IRS wavelength range, we bin 
the data every five datapoints and average their flux values, 
which gives us 30 photometric data points per source. 

We use a x 2 minimisation method to select the opti- 
mum combination of a starburst and an AGN template to 
fit our data. The results can be seen in Fig. [2] We calcu- 
late the 8 — 1000 /im luminosity by integrating the best-fit 
SED; the results are listed in Table[3j We note that six of 
the GOO DS-N sources have the ir infrared luminosities pub- 
lished in |Murphy et al. ( 2009 ) , where a similar procedure 
was used in their calculation. Our L m values agree withi n 
a factor of 2 with the Lm values of Murphy et al. ( 2009 ) . 
Most of the sources of Table[3] require both a star-forming 
and an AGN template to fit the photometric points. The 



1 www.iasf-milano.inaf.it/ - polletta/templates/swire_templates.html 
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spectral decomposition (Fig. [2| suggests that the bulk of 
the mid-IR emission is produced by the AGN component. 
The exception is GSJRS-42. In the case of this source, an 
F-test gives only an 8% probability that the AGN compo- 
nent is needed. 



3.2.4. L x /L 6 ratio 

For faint sources such as those in the CDFs, it is difficult to 
apply X-ray spectroscopy diagnostics owing to the limited 
photon statistics. Alternatively, the Lx to L 6/jm ratio could 
be used to identify any highly obscured, Compton-thick 
AGN. One of the most reliable proxies of the intrinsic power 
of an AGN is considered to be the mid- I R 6 jit m luminosity 
(e.g. [Lutz et al.[ |2004| |Maiolino et~aT| |2007[ ). This wave- 
length region can contain significant emission from the hot 
dust heated by the AGN, and thus provides a re liable diag- 
nostic of the AG N power. On the basis o f this, Alexander] 



et al. (|2008j) and Goulding et al. (20111 claim that very 
low X-ray to L 6 ^ ratios, typically 30-40 times lower than 
those of unobscured AGN, is indicative of very strong at- 
tenuation at X-ray wavelengths and thus of Compton-thick 
sources. 

We present the 2-10 keV luminosity, uncorrected for ab- 
sorption, against the monochromatic 6 /Ltm IR luminosity 
for both the 12 /im sample and the higher redshift sources 
in Fig. [3} We define very roughly a region where Compton- 
thi ck AG N reside, by scaling down the relation of |Fiore et| 
aD([2009|) for the AGN in the COSMOS fields by 0.03. This 



is ver y roughly the am ount of reflected emission in a type-2 
AGN ( |Comastri 2004 1. Moreover, we show the scaled down 
Lx — Lg um relation, by the same amoun t, of AGN in the 
local Universe found by Lutz et al. (2004 1. Finally, we show 
the lu minosity-dependent line predicted by |Maiolino et al.| 
(2007). The luminosity error bars do not play a major role 
relative to the uncertainties introduced by the use of dif- 
ferent Compton-thick lines. We present a typical error-bar 
for one source. This includes both the photon statistics er- 
ror as well as the uncertainty in the spectral model used 
for the conversion from flux to luminosity. Therefore, the 
uncertainties in this diagnostic are quite large. 

Using as a guide our AGN in the local Universe for 
which we have reliable X-ray spectroscopy, we see that 
there is one Compton-thick source (NGC 7582) that lies 
above the local Compton-thick AGN regime. Moving to 
our AGN at higher redshift, there are four GOODS sources 
that could be classified as Compton-thick according to the 
Lx/Le^m diagnostic, i.e. GNJRS-29, GNJRS-55, GSJRS- 
42, and GNJRS-60. Taking into account only the X-ray 
detections, the fraction of possible Compton-thick sources 
would be 66% (4/6). The X-ray upper limits are all situated 
unconvenicntly above the Compton-thick regime. Under the 
extreme assumption that all upper limits are either asso- 
ciated or not with Compton-thick sources, the fraction of 
Compton-thick AGN could vary between 33% and 83 %. 



3.2.5. Contamination by star-forming galaxies. 

Our selection criterion (low EW PAH) is chosen to min- 
imise the number of star-forming galaxies (non-AGN) in 
our sample. The EW upper limits of our FLS sources are 
consistent with a hig h AGN contributio n, close to 100% 
according to Fig. 5 of Sajina et al. (2007). The absence of 



strong PAH features is unlikely to be an artefact of low 
signal-to-noise ratio spectra. There is no trend for a depen- 
dence of the EW on the mid-IR flux (see Table 3). The 
two brightest FLS sources, with X-ray observations avail- 
able, are among the sources with the lowest EW. Another 
possibility could be that the PAHs are possibly suppressed 
in st ar-forming galax ies with large amounts of absorption 
(e.g. [Zakamska] |2010[ ). However, the region responsible for 
the silicate absorption is most probably locate d between 
the AGN torus and the narrow-line-region(e.g. Netzer et 
|al. [ |2007[ ) . The most solid evidence that these silicates are 
mostly located close to the torus is that they are often seen 
incmission in type-1 AGN and in absorption in type-2. In 
contrast, mo st of the PAH em ission should originate from 
larger scales ( Soifer et al. 2002 ) . Furthermore, inspection of 
the strength of the 6.2/zm PAH feature of the rejected r > 1 
sources in Table 5, shows that these display high PAH EW 
(up to EW~ 1.2 ^m) despite there being optically thick at 
9.7/^m. 

However, the most robust indicator of detection of an 
AGN can be provided by the level of its X-ray emis- 
sion. Although a detection at X-ray wavelengths is gen- 
erally sufficient to confirm the presence of an AGN, the 
GOODS fields probe extremely faint fluxes where a num- 
ber of sources could be associated with normal galaxies (i.e 



witho ut an AGN) ( |Georgakakis et al.[ |2007 |Ptak et aL 
2007|. These also have low Lx/L 6/im ratios because their 



X-ray emission is intrinsically weak (instead of absorbed) 
We set the AGN detection limit to Lx — 10 42 ergs -1 as 
no local pure star-forming galaxy has ever presented a 2 



10 keV lu minosity above this threshold (see Tzanavaris et 
all [2006 1. One of the most X-ray luminous star-forming 



ga axies is NGC 3256, which has a 2-10 keV luminosity of 
2.5 x 10 41 erg s -1 , b ut displays no evidence of AGN activity 
(Moran et al. 1999). Therefore, for all our six sources with 



X-ray detections, we can be confident that they are AGN, 
as they present X-ray luminosities above 10 42 ergs -1 . 

In agreement with the above classifications, two of our 
sources (GNJRS-19 and GNJRS-55) have additional evi- 
dence of an AGN via excess of radio emission (Del Moro et 
al., in preparation). FLS- 8196 presents an [Q UI] line with a 
FWHM of w 800 km s" 1 ( [Sajina et aO |2008[ ), wh ich is typi- 
cal of emission from the AGJN narrow-lme-region i Zakamska 



et al.[|2003|). FLS-16080 displays a broad CIVA1550A emis- 



sion line, with a FWHM of ~ 2500 km s , and hence also 



hosts an AGN (Sajina et al. 2008L The presence of a 



broad-line argues against this source being Compton-thick 
There are cases of course of broad-line QSOs ass ociated 
with Compt on-thick AGN but these cases are rare ( |Braito| 



et al. 2004). 



4. Discussion 

4.1. The efficiency of the high-r method for finding 
Compton-thick AGN 

Silicate absorption features in a few (U)LIRGS in the lo- 
cal Universe known t o be associated with Compton-thick 
AGN (e.g. NGC 6240; |Armus et al.] " 
investigate whether a 



20061, prompts us to 
righ optical depth at 9.7 /im could 
provide a reliable diagnostic for the presence of these nu- 
clei. We can reliably identify Compton-thick AGN in the 
nearby Universe by means of X-ray spectroscopy, owing to 
their high X-ray brightness. We argue that six out of nine 
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Table 3. High-r 9.7 /xm AGN at higher redshift 



Name 


a 


S 


z 


s s 




PAH EW 


log[L 6 ] 


log[L m ] 


AGN 


Comment 


AGN 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


CDF-N 


GNJRS-19 


12:35:55.1 


+62:09:01 


1.875 1 


0.097 


2.43 


<0.18 


11.37 


12.62 


0.75 


a,b,d,e 


Y 


GNJRS-29 


12:36:56.5 


+62:19:37 


2.2 2 


0.048 


2.54 


<0.11 


11.83 


12. 62 1 


0.91 


a,d,e 


Y 


GNJRS-30 


12:37:26.5 


+62:20:26 


1.76 2 


0.175 


3.6 


<0.07 


11.92 


12. 74 1 


0.96 


a,d,e 


Y 


GNJRS-55 


12:36:46.7 


+62:14:46 


2.004 3 


0.027 


1.8 


<0.09 


11.49 


12.40 


0.58 


a,b,d,e 


Y 


CDF-S 


GSJRS-14 


03:32:14.5 


-27:52:33 


1.87* 


0.026 


2.22 


<0.20 


11.16 


12.66 1 


0.67 


d,e 


P 


GSJRS-42 


03:32:25.9 


-27:47:51 


1.88 2 


0.02 


1.78 


<0.28 


10.99 


12. 25 1 


0.32 


a.d 


Y 


GSJRS-60 


03:32:40.0 


-27:47:55 


2.0 2 


0.045 


1.35 


<0.17 


11.78 


13.25 


0.57 


a,d,e 


Y 


FLS 


FLS-8196 


17:15:10.2 


+60:09:54 


2.59 2 


12.5 


1.32 


<0.05 


12.50 


13.40 


0.74 


d,e, f 


Y 


FLS-78 


17:15:38.2 


+59:25:40 


2.65 2 


12.7 


1.49 


<0.04 


12.71 


13.51 


1 


d,e 


P 


FLS-8245 


17:15:36.3 


+59:36:14 


2.70 2 


24. 


2.13 


<0.05 


12.38 


13.25 


0.98 


d,e 


P 


FLS-16080 


17:18:44.8 


+60:01:15 


2.01 2 


71. 


2.17 


<0.08 


11.76 


12.85 


0.71 


c,d,e 


Y 


FLS-8550 


17:18:14.6 


+59:56:05 


0.87 2 


255. 


2.74 


<0.27 


11.10 


12. 09 1 


0.70 


d,e 


P 



Th e colum ns are: (T) Name (2), (3 ) Equatorial Spitzer coordin ates (4) Reds hift 



aL] ( |2008[ > or |Murphy et al.| |2009) in the case of GOODS, or |Sajina et~al 



Chapman et al. 



( 2007 in the case of 



(2005) 



; 2 IRS from 


Pope et 


Barger et al. 


(20081. 



case of FLS-8550 where the 11.3 [im EW is quoted). (8) Logarithm of the vL v monochromatic luminosity at 6 fim in units of 
solar luminosity as derived by IRS spectroscopy. (9) Logarithm of the 8-1000 /im IR luminosity in units of solar luminosity as 
derived from the SED fitting. 1 : No photometry at rest-frame wavelengths > 30 /im which constrains the far-IR luminosity. 
(10) Fraction of AGN contribution at 6/im according to the spectral decomposition. (11) Evidence for the presence of an 
AGN: a) X-ray luminosity b) radio-emission c) broad optical line d) absence (or weakness) of PAH feature e) SED f) broad-ish 
(> 500km s _1 ) [OIII] line (12) AGN classification: (Y) secure on the basis of the X-ray luminosity or optical spectroscopy, (p) 
probable on the basis of the mid-IR diagnostics i.e. low-EW PAHs and SED fitting. 



Table 4. X-ray properties of the high-redshift sample 







F,N H 


free 


T = 1.8 


Name 


Lx 


r 


N H 


N H 


(1) 


(2) 


(3) 


(4) 


(5) 


CDF-N 


GNJRS-19 a 


43.5 


1 - Vi -0A4 


i2.i±S:S 


16.8^ 


GNJRS-29 


42.3 






< 7.9 


GNJRS-30 6 


43.1 


54+ 1 - 18 


<17.5 


17.0±1 6 


GNJRS-55 C 


42.3 


2 15+ 26 


<5.8 


<4 


CDFS-S 


GSJRS-14 


<42.9 








GSJRS-42 d 


42.3 






<1.3 


GSJRS-60 


42.1 








FLS 


8196 


<43.8 








78 


< 44.1 








8245 


<43.8 








16080 


<43.7 








8550 


<43.7 









The columns are: (1) Name 
( d ) source 240 in 



Luo et 



( c ) sources 44, 423, and 243 respectively in the catalogue of 



Alexander et al. 



(2003); 



12008). (2) Logarithm of the X-ray Luminosity (2-10 keV) (or 3<r upper limits) uncorrected for 
absorption in units ergs _± . Luminosities have been estimated from the X-ray spectral fits. Where there is no spectral fit, we 

(4) Rest-frame column density in units of 10 22 cm -2 . (5) Rest-frame column density in 



assume T = 1.4. 
units of 10 22 cm" 



(3) Photon index. 
2 , where T = 1.8. 



highly absorbed sources at 9.7^m (rg.7 > 1), in the 12 /im 
IRAS Seyfert sample, with X-ray spectra available in the 
literature are probably associated with Compton-thick nu- 
clei. In any case, all nine sources are heavily obscured AGN. 

Since these are nearby AGN (z < 0.06), they should 
contribute only a small fraction to the X-ray background. 



The peak to the X-ray background is produced at z ~ 
0.7—1 from objects with intrinsic luminosities around 
the 'kne e' of the luminosit y function, i.e., ~ 10 43 — 10 44 
ergs -1 (Gilli et al. |2007 ). The median X-ray luminos- 
ity (uncorrected for absorption) of our six local Compton- 
thick AGN is ~ 0.5 x 10 42 ergs -1 . Assuming a correction 
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Table 5. High-r 9.7 /im sources not included in our sample because of a high PAH EW 



Name 
(1) 



(2) 



(3) 



(4) (5) 



PAH EW log[L 6 ] log[Lx] 



(6) 



(7) 



GNJRS-4 
GNJRS-11 
GNJRS-13 
GNJRS-14 
GNJRS-15 
GNJRS-16 
GNJRS-20 
GNJRS-21 
FLS-283 



12:36:53.4 
12:36:21.3 
12:37:34.5 
12:36:22.5 
12:37:11.4 
12:36:37.0 
12:36:03.2 
12:36:18.3 
17:14:58.3 



+62:11:39 
+62:17:08 
+62:17:23 
+62:15:44 
+62:13:31 
+62:08:52 
+62:11:10 
+62:15:50 
+59:24:11 



1.2T 
1.992 1 
0.64 2 
0.639 2 
1.99 2 
2.02 2 
0.64 2 
2.00 2 
0.94 2 



1.12 
1.55 
1.37 
3.00 
1.09 
3.97 
2.32 
1.96 
1.12 



0.67 
1.18 
0.48 
0.63 
0.78 
0.57 
0.60 
0.66 
0.68 



11.40 
11.41 
11.27 
10.39 
10.74 
11.27 
10.68 
11.19 
11.46 



41.8 
<42.8 
41.3 
41.4 
42.8 
42.1 
41.5 
42.4 
<43.1 



Th e columns are: (1) Name (2) , (3) Equatorial Spitzer coo rdinates (4) Redshif t: Chapman et al. (20051; IRS from Pope et 
aLl ( |2008| ), |Murphy et al.| ( |2009| ) in the case of GOODS or |Sajina et al.| ( |2007|) in the case of PEST (5} optical depth at 9.7/nm. 
TJJ) PAH equivalent-width at 6.2/tm, in units of fim (except in the case of FLS-283 where the 11.3 /im EW is quoted). (7) 
Logarithm of the vL v monochromatic luminosity at 6 fim in units of solar luminosity as derived by IRS spectroscopy. (8) X-ray 
luminosity (uncorrected for obscuration) in the 2-10 keV band in units ergs -1 . 



of ab out 33 for the intrinsic X-ray luminosity (Comastri 
20041 , this translates to « 2 x 10 43 ergs -1 . Hence, AGN 
with comparable luminosity, but at higher redshift, are the 
typical AGN responsible for the peak of the X-ray back- 
ground. 

At higher redshift, the fraction of candidate Compton- 
thick AGN is more uncertain. This fraction can be derived 
indirectly from the X-ray to IR luminosity ratio. If we take 
only the X-ray detections into account, this is close to 66%. 
When we also consider the upper limits, we find that the 
uncertainty becomes large: the fraction of Compton-thick 
sources could vary between 33% and 83%. 

The frequency of Compton-thick AGN measured above 
is comparable to or higher than derived by other mid-IR 
methods for finding Compton-thick sources. For example, 
on the basis of X-ray stacking analyses , several authors (e.g. 
Daddi et all |2007[ [Fiore et al.| |2008l |Georgantopoulos et 



al. 2008 1 proposed that the 24 /zm-excess sources are good 
candidates for hosting Compton-thick nuclei. However, be- 
cause of the lack of X-ray spectroscopy the fraction of 
Compton-thick sources cannot be readily quantified. Based 
on brighter sources with X-ray s pectroscop y in the GOODS 



dGeorgan topoulos et al.[ 2011 
al.| |20l0[) and SWIRE fieWj 



AEGIS (jGe orgakakis et 
Lanzuisi et al.[ 2009p 



lEe 

fractions of Compton-thick sources found among IR-excess 
AGN range between 0% and 50%. Another widely applied 
mid-IR method for the selection of Compton-thick AGN 
( [Alexander et al." 2008) relies on low X-ray to IR ratios. 
The principle behind this method is tha t a low Lx/Lir ratio 
is a s i gn of heavy obscuration (but see Yaqoob fc Murphy] 
20101 iGeorgakakis et all 120101). 



Although we have shown that the success rate of our 
technique in finding Compton-thick AGN is close to 70%, 
at least among optically confirmed Seyferts in the local 
Universe, there are two caveats. First, the detection of this 
feature alone does not imply the presence of an AGN in 
these systems. There are star-forming systems with large 
absorption features at 9.7/xm that are associated instead 
with star-fo rming galaxies. One case is NGC 3628 for which 
Tg.g = 1.64 ( |Brandl et al. 2006), although the X-r ay emis- 
sion in this system does not come from an AGN (Dahlem 
et al.| |1996[ ). In other words, the high 79.7 optical depth 
characterises systems with large amounts of dust but does 



not necessarily prove the presence of an AGN. The sec- 
ond caveat has to do with the completeness of our method. 
There are well-known Compton-thick sources in the litera- 
ture that contain no significant 9.7 /im absorption. We dis- 
cuss this in detail in the following section. 



4.2. Si optical depth versus X-ray column density 



Shi et al. (2006) described observations of 9.7 /im features 
in 97 AGN including a wide range of types. They found that 
the strength of the silicate feature correlates with the hy- 
drogen column density, derived from X-ray observations in 
the sense that the low column densities correspond to sili- 
cate emission, while high columns correspond to silicate ab- 
sorption. They point out that the column densities derived 
from X-ray spectral fitting are always higher than those 
estimated on the basis of the Si optical depths. Moreover, 
there appears to be a large scatter between the X-ray and 
the 9.7 /im column density in the sense that the same 79.7 
corresponds to X-ray column densities that could differ by 
as much as an order of magnitude. This result can be seen 
more clearly in Fig.|4j where we plot the 9.7 /im optical 
depth as a function of the X-ray column density for bot h 
the 12 /jm sample as well as the sample of Shi et al. (2006). 
We note that the latter is based mainly on th e optic ally se- 
lected sample of Risaliti, Maiolino & Salvati (19991. Some 
column d ensities in this sample were revised by |Malizia et| 
[aL] ([2009 ) and we include these updated values in Fig.UJ 
This figure clearly shows that the majority of the optically 
thick sources detected in the mid-IR display Compton-thick 
absorption at X-ray wavelengths. In addition we see that 
there are numerous Compton-thick AGN according to X- 
ray spectroscopy that have optical depths r g 7 < 1. 

Perhaps the most well studied case of a low Si r 
Compton-thick source is NGC 1068. The X-ray column den- 
sity of this Seyfert-2 galaxy is Nh > 10 25 cm" 2 as derived 



from high energy X-ray spectroscopy ( Matt et al. 2004 ) 



with the BeppoSAX mission. However, the mid-IR spec 



trum shows no prominent 9.7 /im absorption (Sturm et al. 



2000). In Shi et al. (2006) there are many more examples 



of Compton-thick nuclei which are optically thin at 9.7 /im 
(see Fig. |4]). From the above discussion, it becomes clear 
that although the presence of a high— r absorption feature 
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Fig. 2. Mid-IR rest-frame SEDs of the GOODS and FLS sample. The blue and red line denote the torus and star-forming 
contribution while the purple line denotes the sum of the two. 



implies the presence of a Compton-thick AGN, with a prob- 
ability of roughly 70%, the opposite is not true, i.e., there 
are many Compton-thick AGN which are not significantly 
absorbed at 9.7 /j,m. 

The reason for this could be that the X-ray and in- 
frared emitting regions are not physically coincident and 
the infrared arises from a region further away from the ac- 



cretion disk. It is possible that within the framework of a 



clumpy torus model (see e.g. Nenkova et al. 2008[ Wada, 
Papadopoulos fc Spaans] I2009| |Nikutta et al.[|2009p , the" 
X-ray emission from the AGN could be obscured by a cir- 
cumnuclear cloud, while the infrared emission, which is 
produced at larger distances from the AGN, c ould expe 
rience low or no obscuration at all (see Fig. 5 in Shi et al. 
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Fig. 3. Rest-frame (uncorrected for absorption) Lx versus Lg^m luminosity diagram. The filled triangles correspond to 
the Compton-thick AGN in the f2/im sample, while the open triangles correspond to the three heavily obscured but 
not Compton-thick AGN in the same sample. The red (blue) filled circles or upper limits denote the high-r sources 
in the CDF-N (CDF-S). The black upper limits denote the FLS points. Finally, the green points (filled squares and 
small upper limits) denote the sources with considerable PAH emission that have not been included in our sample. The 
dash line (magenta) denotes the low Lx/L6 Mm area that should be popula ted by Compton-th ick AGN. This line is 
derived from the average Lx/L 6Aim relation derived for the COSMOS AGN (Fiore et al. 2009) by scaling it down by 
a factor 0.03, which corresponds to the average flux of the reflection compo nent in the 2-10 keV band relative to the 



intrinsic power-law component in local Compton-thick AGN (Comastri 2004 ). The solid line (green) corresponds to the 
luminosity-dependent Compton-thick line based on thelMaiolino et al. (2007) relation, scaled down by the same amount. 



Finally, th e horizontal dotted line at Lx 
processes (Tzanavaris et al. 2006). 



10 ergs 1 denotes the upper limit to X-ray emission due to star-forming 



2006). |Spoon et al.|p007| proposed that the depth of the Si 
9. 7/Ltm feature can be used as a tool to predict the obscur- 
ing dust distribution. In this scenario, high silicate depths 
correspond to sources embedded in a smooth dust distri- 
bution, while low depths correspond to those covered by 
patchy absorbers. 



4.2.1. Comparison with previous results 

We to compare our results with the work of Alexander et 



al. (2008), who searched for Compton-thick AGN in the 



CDF-N using both IRS spectroscopy in com bination with 



the Lx/L 6Mm ratio. Alexander et al. (|2008) discussed the 



properties of two sources with available mid-IR IRS spec 
troscopy in the CDF-N, claiming that both could be as- 
sociated with Compton-thick AGN. One (HDF-oMD49) is 
an optically identified AGN that is formally undetected in 



the 2 Ms data [Alexander et al.| ( |2003[ ). The Spitzer IRS 
spectrum and spectral energy distribution of this object is 
AGN dominated. On the basis of the source's low Lx/Lir 

(120081) 



ratio, 



Alexander et al. 



argue that this AGN is 
Compton-thick. The otherliource (SMMJ123600+621047) 
is a z=2.002 submilli meter-emitting gal axy (SMG) with a 
mid-IR-bright AGN ( |Pope et al~| [2008) that is again un- 
detected in the 2 Ms CDF-N data. The Lx/L 6/im luminos- 
ity ratio indicates that it hosts a Compton-thick AGN, al- 
though optical and near-IR spectroscopic observations do 
not reveal the signatures of AGN activity, leaving open the 
possibility that this source is associated with a star-forming 
galaxy and not with an AGN. None of these sources display 
a t > 1 9.7/zm absorption feature in its IRS spectrum. 
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tcctcd in the X-rays, while all five FLS sources remain 
undetected. For the X-ray detections, the high X-ray 
luminosities (> 2 x 10 42 ergs _1 ) confirm that all these 
sources are AGN. For two FLS sources, their optical 
line emission also indicates that an AGN is present. 
Because our sources are faint in X-rays, we are un- 
able to ascertain from X-ray spectroscopy whether at 
least some of our sources are Compton-thick. The use 
of the Lx/L 6Aim ratio indicates that the success rate of 
the high-r method in identifying Compton-thick AGN 
could be about 66 % (4/6), if we consider X-ray detec- 
tions only. However, it can be as low as ~33 % (4/12), 
or as high as 83% under the extreme assumptions that 
all the X-ray upper limits are either associated or not 
associated with Compton-thick sources. Deeper X-ray 
spectroscopy will be needed to properly constrain the 
exact fraction of Compton-thick AGN among high-z, 
high-r sources. 

In the future, observations with the JWST will be able to 
10 22 1 23 1 24 1 25 1 26 identify numerous galaxies with deep Si features that could 
^x-rays ( cm -a) be used to compile large candidate Compton-thick AGN 

samples. 




Fig. 4. Optical dept h at 9.7 fj,m vers us X-ray hydrogen col- 
umn density for the |Shi e"t~ al. (2006) sample (black points) 
and the optically-thick AGN in the IRAS sample of Wu 



et al. (2009) (red points). The solid line marks the region 
where the X-ray Nr equals that inferred from T9.7, assum- 
ing a Milky- Way extinction curve. 



5. Summary 

The goal of this work has been to investigate whether the 
presence of high-absorption at mid-IR wavelengths (T9.7 > 
1) can be used as an efficient means for identifying highly 
obscured or even Compton-thick AGN. Our conclusions can 
be summarised as follows: 

— Our re-analysis of Spitzer-IRS spectra shows that there 
are eleven AGN with rg.7^ m > 1 in the 12/im Seyfert 
IRAS sample. For nine sources X-ray spectroscopic ob- 
servations are available in the literature. The X-ray 
spectra show that a large fraction (six of nine) is prob- 
ably associated with Compton-thick AGN. The success 
rate of this method is comparable to or even higher than 
those of other methods using either mid-IR photometric 
information or the X-ray to IR luminosity ratio. 

— Even though optical thickness in the MIR, as deter- 
mined by 9.7 micron Si absorption with r > 1 can be 
efficiently used to identify Compton-thick sources in the 
local Universe with a high success rate, this method can- 
not provide complete samples of these sources. This is 
because a large fraction of Compton-thick sources do 
not have large high Si optical depths. This could pos- 
sibly be explained within the framework of a clumpy 
torus model. 

— Using Spitzer-IRS spectra in the CDF and the FLS sur- 
veys, we have compiled a sample of 12 sources in total, 
with r 9 .7 > 1 at higher redshift (z=0.87-2.70). These 
are selected to have weak PAH features to maximise 
the probability that they host an AGN. The vast ma- 
jority of sources in the CDFs (six out of seven) are de- 
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